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Abstract 

The recently published experimental data for specific heat C p of liquid 
helium in zero gravity conditions very close to the A-transition have been 
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J> . discussed. We have shown that these data allow different interpretations. 

^£ | They can be well interpreted within the perturbative RG approach and within 

our recently developed theory, as well. Allowing the logarithmic correction, 
the corresponding fits lie almost on top of each other over the whole range 
of the reduced temperatures t (for bin averaged data) 6.3 • 1CT 10 < t < 
q | 8.8 • 10~ 3 . However, the plot of the effective exponent a e s(t) suggests that the 

behaviour of C p , probably, changes very close to T\. To clarify this question, 
we need more accurate data for t < 10~ 7 . In addition, we show that the 
experimental data for superfluid fraction of liquid helium close to T\ within 
t € [3 ■ 1CT 7 ; 1CT 4 ] can be better fit by our exponents v = 9/13 ~ 0.6923, 
A = 5/13 ~ 0.3846 than by the RG exponents v ~ 0.6705 and A ~ 0.5. The 
latter ones are preferable to fit the whole measured range t G [3 • 10~ 7 ; 10~ 2 ] 
where, however, remarkable systematic deviations appear. Our estimated 
value 0.694 ± 0.017 of the asymptotic exponent v well agrees with the theo- 
ry \ retical prediction v = 9/13. 
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1 Introduction 

It is widely accepted to consider the measurements in liquid helium near A-transition 
point T — T\ as a crucial test of validity of the theoretical predictions for the critical 
exponents, since these measurements are done with a high degree of accuracy much 
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closer to the critical point than in any other experiments or numerical simulations. 
In particular, it is believed that accurate experimental measurements of specific heat 
C p of liquid helium very close to the A-transition point in zero gravity conditions (in 
space) PP provide a convincing evidence of overall correctness of the perturbative RG 
approach. The aim of our paper is to show that this conclusion is not unambiguous, 
since these experimental data as well as those of the superfluid fraction of liquid 
helium can be equally well or even better interpreted by a completely different set 
of critical exponents provided by our recently developed theory 

2 Interpretation of the specific heat data 

It has been found in pQ that fits of experimental data for a wide range of reduced 
temperatures 5 ■ 1CT 10 <t< 10~ 2 below T A by using two slightly different ansatz, 

C p = ^-t~ a (l + a-t A + b-t 2A ) + B~ (1) 

and 

C p = ^r a (l + a~t A ) + b~t + B~ , (2) 

biased by the RG theoretical value of the correction-to-scaling exponent A = 0.529, 
provide well consistent values of the specific heat exponent a = — 0.0127 ± 0.0003 in 
a good agreement with the value —0.01294 ± 0.0006 of the variational perturbative 
theory |3] as well as in a worse, but still acceptable, agreement with more recent 
estimates a = -0.01126 ± 0.0010 | and a = -0.0146 ± 0.0008 0. Apart from 
the exponent a, some other quantities have been determined and compared with 
the RG values in PQ. However, the agreement is not so good to conclude that any 
theoretical approach, which does not agree with the perturbative RG, is wrong. In 
particular, the experimental quantity P = (1 — A + /A~)/a is 4.154 ±0.022, whereas 
the recent RG calculation (Ref. 63 in [T]) yields P = 4.433 ± 0.077. 

We have found that the measured data of C p PQ can be well reproduced also by 
an ansatz of the form 

C p = t- a {C + Alnt)(l + at A )+B, (3) 

with fixed exponents a = —1/13 and A = 5/13 proposed in [21 El- It is consistent 
with the idea that specific heat can have a logarithmic correction, as discussed in j2] . 
The power-like singularity is recovered at A = 0. Note that in [2] (cf. Eq. (60) there) 
an ordinary term ~ t~ a is related the behaviour of the correlation function within 
the range of wave vectors k ~ l/£, where £ ~ t~ u is the correlation length, whereas 
the logarithmic term can appear due to the contribution of the region k ^> l/£. In 
this aspect, the ratio A/C in Q can be varied in a wide range of values. 

From the raw data of pQ given in [7] we have produced the set of bin averaged 
data points by dividing each decade of the reduced temperature t = 1 — T jT\ in 
10 segments of equal width when looking in the logarithmic scale. One binned 
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data point has been obtained by an averaging over C p and t values within one 
segment, and only the data for the smallest t values within a twice wider interval 
4.7 ■ 1CT 10 < t < 7.9 ■ 10 -10 have been merged together into one bin to reduce the 
statistical error. In our binning, the averaged data points come as close to T\ as 
t > 6.3 ■ 1CT 10 , whereas those given in [2j extend only to t > 7.94 ■ 1CT 10 . 

The percent deviations from the least-squares fits to (jUJ) and (JJJ) are shown in 
Fig. [TJ The upper picture corresponds to © with fixed exponents a = —1/13 and 
A = 5/13 and coefficients C = -167.536, A = 11.6593, a = 0.19788, and B = 
198.26, whereas the lower one represents the fit to (0) with exponents a = —0.01264 
(fit parameter) and A = 0.529 (fixed) and coefficients listed in Tab. II of pQ. As 
in some fits made in pQ, we have assigned the error bars 0.02% to the bin averages 
which originally had smaller errors. In this way, we have reduced the impact of these 
data points, located at relatively large values of the reduced temperature t, where 
the asymptotic ansatz © is not very accurate. As we see from Fig. [JJ both fits are 
almost identical in the whole range of the reduced temperatures. The fit with our 
exponents (top) is slightly worse at the largest t values. It can be well understood, 
since (to ensure the stability of the fit parameters) we have neglected the subleading 
correction of the kind t 2A included in the other ansatz (0). Besides, our fit is even 
slightly better at the smallest t values: the mean percent deviation for 10 smallest 
t values is -0.425 ± 0.690 in our case of © and -0.975 ± 0.686 in the case of ((TJ). 
These deviations are reduced to 0.004±0.695 and — 0.477±0.691, respectively, when 
shifting the T\ — T values by 0.5 nK within the experimental error bars pQ. 

Note that only the possibility and not the necessity of the logarithmic correction 
follows from the theory [2] . However, the presence of the logarithmic correction for 
specific heat, perhaps, is a quite general feature: the logarithmic singularity (as a 
special case of the logarithmic correction when a = 0) of specific heat is a rigorously 
stated fact in 2D Ising model and our Monte Carlo simulation data for 3D Ising 
model jH] also supports the logarithmic singularity. Our analysis of the experimen- 
tal data for the superfluid fraction, made in Sec. H3 gives one more argument: it 
suggests that the exponent v is remarkably larger than 0.6705. According to the 
known scaling relation a + dv = 2 jS], the exponent a then should be remarkably 
more negative than —0.0115. It well coincides with the measured data £[| in the 
whole range of the reduced temperatures only if the pure power is perturbed by a 
logarithmic correction, as proposed by ansatz (jSJ). 

Considering a as a fit parameter in (JHJ), we obtain a value a = —0.0848 ±0.0039, 
which is quite close to our theoretical prediction a = —1/13 ~ —0.0769. The small 
systematic deviation could be caused by the error of the asymptotic ansatz Q at the 
largest t values. This problem cannot be reliably solved by adding more correction 
terms or by narrowing the range of the fit, since the minimum of the \ 2 f° r such fits 
is very shallow, i. e., the results become poorly defined. 

Alternatively, we have fit the data within a moving window t G [tf, lOOtJ of the 
reduced temperatures to the simplest possible ansatz 

C p = At- a + B, (4) 
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Figure 1: Percent deviation of the fitted C p data points from the ansatz (jBJ) with 
fixed exponents a = —1/13 and A = 5/13 (top) and from the ansatz (JIJ with the 
fit parameters found in (bottom). 
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Figure 2: The effective exponent o^g vs t A estimated by fitting the bin averaged 
data within [tf, lOOtj] to the ansatz (£Q) with unbiased (left) and shifted according 
to T\ — > T\ + 0.5 nK (right) values of the reduced temperatures. The horizontal 
dot-dashed line indicates the a value —0.01264 found in jT], whereas the dashed 
curve shows schematically a behaviour expected from [2]. 

where t« corresponds to the i-th bin averaged data point. It yields the effective 
exponent a^(t), where t belongs to the considered interval. For convenience, we 
have defined it as t = y/t m i n t max , where t min = £j is the minimal and t max is the 
maximal t value in the interval. By this method, the result converges to the true 
value of a at t — > irrespective to the error of the asymptotic ansatz at finite 
t. It works also when the logarithmic correction is present, only in this case the 
convergence is very slow, like a — a e s{t) ~ 1/lnt at t — *■ 0. The only problem is 
that this method requires that the measurement errors both for C p and t remain 
sufficiently small when approaching T\. 

We have plotted in Fig.|2]the results for a eff depending on t A (with our exponent 
A = 5/13) obtained by using the unbiased values of (left), as well as the shifted 
values t' i = t i + 2.3 • 10~ 10 (for the same sets of measured points) corresponding to 
T\ — > T\ + 0.5 nK (right). In spite of the large error bars, these plots show certain 
trend, where the effective exponent a e s(t) tends to decrease below the (RG) value 
—0.01264 found in [Tj. Moreover, we have verified that the same trend is observed 
in both cases when only the odd and only the even raw measurements (the original 
values listed in [7|) are used. Hence, we cannot exclude any striking scenario. For 
instance, a e s(t) could converge to our asymptotic value a = —1/13, as indicated by 
dashed lines, particularly, if we allow a small shift in T\ — T values (right picture) 
within the experimental error bars. However, such a behaviour would mean that 
the logarithmic correction is absent, since the convergence is rather fast. It would 
imply also that the estimation of a from the fit over the whole measured range is 
not valid: formally, such a fit looks good, but it effectively ignores the systematic 
deviations at the smallest t values where the error bars are larger. 
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Due to the experimental errors, the results of our analysis of the effective expo- 
nent are not conclusive, they only point to a possible scenario. From an intuitive 
point of view, it does not seem plausible that such a remarkable change in the be- 
haviour of the system could take place at so small reduced temteratures (t < 1CT 7 ). 
However, since the deviations from the asymptotic scaling law are caused by the 
corrections to scaling, an essential parameter is t A rather than t, and the values 
of t A in Fig. H are not so extremely small. Besides, the critical region, where an 
asymptotic ansatz is valid, can be as narrow as t A ~ 0.001 even in a simple mean 
field model. An example is given in 9J (p. 75). From this point of view, it is possible 
that the deviations in Fig. El represent a real physical effect and not an artifact. On 
the other hand, random deviations in Fig. too often are as large as 2 standard 
deviations or even larger, therefore the unusual behaviour of a e s in Fig. El can be 
ascribed also as an artifact. 

3 Interpretation of the experimental data for 
the superfluid fraction 

Here we discuss the experimental data for the superfluid fraction p in liquid He. It 
decreases asymptotically (at t — > 0) as p ~ tt. It is believed (see (TUj and references 
therein) that the exponent ( is equal to the correlation length exponent v for the 3D 
XY model. In j2] , the superfluid fraction of 4 He measured in [TU] has been discussed 
with an aim to compare the experimentally observed behaviour at the temperatures 
closest to T\ with our theoretical prediction v = 9/13 

The data listed in [TTj allow a more precise comparison. For this purpose, first 
we have fit these data to the asymptotic ansatz 



including two corrections to scaling. Similar fits over the whole measured range 
t G [3 ■ 10~ 7 ; 10~ 2 ] have been considered in (TUl [Uj. Note that at A = 0.5, used 
in [TUJ HI], the second order correction reduces to the analytical one, and (jSJ) differs 
from the ansatz of [TO*] ITT] only by a remainder term of higher order. The overall 
fits discussed in |10[ ITT] yield v ~ 0.6705 in agreement with the RG prediction and 
in disagreement with our value v = 9/13. However, these fits look really good only 
within t G [10 -5 ; 10 -2 ], whereas remarkable systematic deviations appear at smaller 
t values. This phenomenon was discussed in ^JJ and no reasonable explanation was 
found. In particular, the effect of gravity is negligible in these experiments jTTJ and 
the ±20 nK uncertainty in T\ — T also does not explain these systematic deviations. 

Our theory [2] provides an explanation. First, the data cannot be well fit within 
the whole measured range t G [3 ■ 10~ 7 ; 10" 2 ] simply because the critical region, 
where the asymptotic expansion in t powers is valid, is much narrower than 10~ 2 . 
Second, the data can be well fit to (0) with our exponents v = 9/13 and A = 5/13 
within a reduced range t G [3 • 10~ 7 ; 10~ 4 ], which means that the measured data for 
t < 10~ 5 are not anomalous, but the region of validity of © is as narrow as 10 -4 . 




(5) 
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Figure 3: Percent deviation of the experimental p (superfluid fraction) data points 
from the least squares fit to ansatz (j5J) within t E [3 • 10 -7 ; 10~ 4 ] with our exponents 
v = 9/13, A = 5/13 (left) and with the RG exponents v = 0.6705, A = 0.5 (right). 

The percent deviations from the least-squares fits within t 6 [3 • 10 -7 ; 10" 4 ] with our 
(left) and RG (right) exponents are shown in Fig. El As we see, in our case there are 
no essential systematic deviations, whereas in the RG case they are observed like in 
the case of the fit over the whole measured range [HH ITTj . 

Similarly as in Sec. |2l we have evaluated also the effective exponent u^(t) as the 
local slope of the hip vs hit plot within [tf, 5U], where U is the reduced temperature 
of the i-th measurement and t is the middle point of the fitted interval in the 
logarithmic scale. The results depending on t A (with our value A = 5/13) are 
shown in Fig. 0] Evidently, the effective exponent tends to deviate above the value 

0. 6705 (dot-dot-dashed line) obtained in [T0J[TT]. On the other hand, the fit of this 
plot to a parabola (solid line) gives the asymptotic estimate v = 0.694 ± 0.004 in 
excellent agreement with our theoretical value 9/13 ~ 0.6923. However, taking into 
account the ±20 nK uncertainty in the T\ value JU], the error bars become larger, 

1. e., v = 0.694 ±0.017. 

4 Conclusions 

Although the opinion dominates in publications that the perturbative RG theory 
is strongly confirmed by very accurate measurements of the specific heat and the 
superfluid fraction in liquid helium near the A-transition point, our current analysis 
shows that these experimental data can be well understood and interpreted also 
within our recently developed theory [2|. 
In summary we conclude the following: 

1. The critical exponents of the perturbative RG theory look preferable from a 
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Figure 4: The effective exponent u e s vs t A determined from the local slopes of the In p 
vs hit plot. The lower dot-dot-dashed line indicates the (RG) value 0.6705 obtained 
in |10[ lllj. whereas the upper dashe line shows our theoretical value v = 9/13. The 
solid curve represents the least-squares fit to a parabola. 

point of view that all measured data points must be necessarily fit on one 
curve and no logarithmic corrections are normally expected. However, if we 
allow the logarithmic correction to specific heat C p , our theory also provides a 
good fit of C p data for the whole measured range. Our fit then is slightly worse 
at the largest reduced temperatures t ~ 10~ 2 and better at the smallest ones 
t ~ 10~ 9 . As discussed in Sec.|21 the existence of the logarithmic correction is 
partly supported by our estimation of the exponent v (Sec. EJ), as well as by 
some general argument. 

2. The analysis of the effective exponent a e fr indicates that the behaviour of C p 
could be remarkably changed very close to T\, in such a way that the true 
asymptotic singularity is power-like (without the logarithmic correction) with 
the exponent which more probably is closer to our value a = —1/13 than to 
the (RG) value —0.01264 obtained in pQ. However, this effect can be ascribed 
also as an artifact caused by the measurement errors. Further improvement 
of the experimental accuracy for t < 10~ 7 would be very helpful to clarify this 
question. 

3. As compared to the RG exponents, our critical exponents are better consistent 
with the closest to T\ data (t G [3 ■ 10~ 7 ; 10~ 4 ]) for the superfluid fraction 
(cf. Figs. Eland |3J). A self consistent estimation in this case yields v = 0.694 ± 
0.017 in agreement with our theoretical prediction v = 9/13 ~ 0.6923. Since 
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the fit over the whole measured range (t e [3 • 10 7 ; 10 2 ]) in no case is really 
good, we argue that our way of estimation is preferable. 
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